Kinks and waterfalls as signatures of competing order 
in angle-resolved photoemission spectra of La 2 _ x Sr x Cu04 



(N 

o 

3 



G. Mazza 1 ' 2 , M. Grilli 2 - 3 , C. Di Castro 2 - 3 , and S. Caprara 2 ' 3 
1 International School for Advanced Studies (SISSA), Via Bonomea 265 34136 Trieste, Italy 
1 CNISM and Dipartimento di Fisica, Universita di Roma "La Sapienza" Piazzale Aldo Moro 5, 1-00185 Roma, Itah 
3 Consiglio Nazionale delle Ricerche, Istituto dei Sistemi Complessi, via dei Taurini, 1-00185 Roma, Italy 

We show that the so-called kinks and waterfalls observed in angle-resolved photoemission spectra 
of La2- a; Sr :E Cu04, a prototypical high-T c superconducting cuprate, result from the coupling of 
quasiparticles with two distinct nearly critical collective modes with finite characteristic wave vectors, 
typical of charge and spin fluctuations near a stripe instability. Both phonon-like charge and spin 
collective modes are needed to account for the kinked quasiparticle dispersions. This clarifies the 
long-standing question whether kinks are due to phonons or spin waves and the nature of the bosonic 
mediators of the electron-electron effective interaction in L^-iSr^CuCU. 
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The metallic phase of high-T c superconducting 
cuprates evolves remarkably with changing the tempera- 
ture T and the doping x. A normal Fermi-liquid behavior 
is only found in overdoped samples, with x larger than 
the optimal value x op t (where the maximum supercon- 
ducting critical temperature T c is achieved) . At x ~ iTopt 
and T > T c the metallic phase seems to be ruled by the 
temperature as the only relevant energy scale, a typical 
signature of quantum criticality. In underdoped sam- 
ples, with x < x opt , an even stronger anomaly is found, 
with a pseudogap opening around the Fermi energy, be- 
low a doping dependent temperature T*(x). Whether 
this is accompanied by the onset of some sort of order- 
ing is still matter of debate. Nonetheless, models with 
nearly critical collective modes (CMs) coupled to fermion 
quasiparticles (QPs) may not only explain the anomalous 
metallic phase, but also provide candidate mediators of 
a retarded pairing interaction (the so-called glue) [IHij], 
alternative to phonons in ordinary superconductors, and 
are therefore actively investigated. Various proposals for 
sources of nearly critical CMs include the antiferromag- 
netic phase at x ~ 0] , time- reversal-breaking plaquette 
currents [j|, order parameters with exotic wave symme- 
try d, 0| , or stripe ordering [1, Q . 

In this Letter, we show that the so-called kinks and 
waterfalls [l(| observed in angle-resolved photoemission 
spectroscopy (ARPES) identify charge (C) and spin (S) 
CMs on the verge of a stripe instability as the main source 
of scattering in La 2 _ x Sr :c Cu04 (LSCO), solving the long- 
standing phonon-vs-spin issue, at least in LSCO. The 
sudden changes in the QP velocity (kinks) occurring at 
different energies in different regions of the Brillouin Zone 
(BZ), have often been attributed to a phonon with a near 
frequency 11 , [l2| , and it has been argued that S fluc- 
tuations could account for the kinks in YBa 2 Cu307_z 
(YBCO) Q, and Bi 2 Sr 2 CaCu 2 8 (BSCCO) Q. We 
show that the same CMs also account for the sud- 
den, nearly vertical, drops of the QP dispersions at 
high/moderate binding energies (waterfalls) 



0.15 — 0.26 and various T provided evidence for two dis- 
tinct CMs peaked at finite characteristic wave vectors in 
LSCO 18]. One CM, essentially propagating and cen- 
tered at typical phonon frequencies, is associated with C 
fluctuations (strongly mixed with the lattice degrees of 
freedom) near an incommensurate charge-density wave 
instability. The other CM, more diffusive and extend- 
ing to higher energies, is associated with S fluctuations 
peaked near the wave vector of antiferromagnetic or- 
der. The behavior of the characteristic low-energy scale 
of the two CMs suggests that a quantum critical point 



occurs at xqcp ~ 0.19 [2lM23| . associated to a phase 



with stripe-like C and S modulation, whose onset occurs 
via a harmonic incommensurate charge-density wave at 
T rj T*(x) @, 0, Remarkably, the strengths of the 
two CMs have an opposite x dependence (see Fig. 5 in 
Ref. [13]): The strength of the S CM decreases with in- 
creasing x and almost vanishes in the most overdoped 
sample (x = 0.26), whereas the strength of the C CM 
increases with increasing x and tends to saturate in the 
overdoped regime. The value x « 0.19 marks the bound- 
ary between the spin- and the charge-dominated region. 
At small x, S fluctuations are naturally enhanced by in- 
cipient antiferromagnetism, whereas in the optimally and 
overdoped regime C fluctuations dominate. 

Here, we phenomenologically proceed to analyze the 
implications of the same two CMs, as derived from Ra- 
man experiments, on QP spectra. We consider the gen- 
eral Gaussian form of CM propagator 

2? A (q,w n ) = -— - 1 , - - - , (1) 
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A previous survey of Raman spectra jl7| for x 



where A = C, S, uj n is the bosonic Matsubara frequency, 
T A (q) = ra\ + v\ [2 — cos(q x a) — cos(q y a)] describes the 
dispersion of a lattice periodic CM, and reproduces the 
behavior T A (q— Q A w0)w to a + \v\(a — Qx) 2 obtained 
in different contextes for C [9| and S [j, [26( CMs. m\ 
is proportional to the inverse squared correlation length 
£^ 2 , v\ sets the curvature at the bottom of the CM dis- 
persion law, and a is the spacing of the two-dimensional 
square lattice describing the Cu0 2 planes of cuprates, 
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henceforth taken as unit length. The propagator ([I) is 
peaked at a characteristic wave vector Qa, has a diffusive 
character at low energy, and becomes more propagating 
above the energy scale The CM dispersion is lim- 
ited by an energy cutoff A\, setting a momentum cut- 
off |q A | w (Aa/^a) 1 ^ 2 - The values of the characteristic 
wave vectors, Qc and Qs, are extracted from neutron 
scattering experiments: The incommensurability of the 
S density modulation, half of the incommensurability of 
the C density modulation when observed at x = 1/8 
(30l l3lj |. saturates for x > 1/8, yielding [32| Qc = 
7r(±i,0),7r(0,±i), Qs=7r(l±i,l),7r(l,l±|). 

We adopt for the fermion QPs on the CuC>2 planes of 
LSCO a tight-binding dispersion law including nearest 
(t = 400 meV) and next-to-nearest (f = — 0.2W) neigh- 
bor hopping terms, 

e k = —It (cos k x + cos k y ) — At' cos k x cos k y — /1, (2) 

where /i is the chemical potential. Similarly to the 
electron-phonon coupling, QPs are here coupled to CMs 
through dimensional coupling constants g\. The survey 
of Raman data on LSCO 17[ yielded the doping evolu- 
tions of m, A, fl, and n = g 2 jtv, reported in Tab. HI Our 
aim is to fit ARPES data with the same CM parameters, 
although it should be borne in mind that Raman response 
is a momentum integrated quantity, so that the precise 
q dependence of the CM-mediated effective interaction, 
and therefore the value of q, is not fully constrained. We 
adjust q as a fitting parameter [27j ]. which in turn fixes 
v = A/q 2 and g = \Jntv. The obtained values (see Tab. 
I), yield four peaks in the CM dispersion, consistent with 
the observation of four separate peaks for the S CM in 
neutron scattering experiments 

The effect of CMs on QP spectra is captured comput- 
ing the lowest-order QP self-energy S(k, uS) = Ec(k, w) + 
Eg(k,a;), and the single-particle spectral density 



|ImS(k,w) 



7T [u - e k - ReS(k, u)f + [ImS(k, lu)} 2 
The imaginary part of the self-energy is 
f d 2 

ImE a (k, a) = g\ / —-^ 1qM~Q>) 
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where a sum over the star of equivalent wave vectors Qa 
is understood, /±(w) = 1/ (e"/ T ± l), and the smooth 
cutoff function 7g(q) = exp [— (2 — cos^ — cosq y ) /Q 2 ] 
accounts for the suppression of the CM-QP coupling away 
from Qa- As an order-of-magnitude estimate, Qa ~ |<5aI- 
For any given k and w, we numerically integrate Eq. 
© and obtain ReE via Kramers-Kronig transformation. 
The chemical potential fi is fixed imposing that 

<i 2 k f +co 

dio A(k, ui)f+(uj) = 1 — x. 
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TABLE I. Parameters of the C and S CMs: m, A, fi, and k 
are extracted from Raman data LL31. 



We calculate the ARPES intensity convoluting A(k, oS) 
with a Gaussian of width w lOmeV, mimicking energy 
resolution, and considering only the occupied states. 
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FIG. 1. (Color online) Spectral density along the cuts of the 
BZ reported in Fig. for x = 0.15 and T = 40 K. Dots 
represent the maxima of the MDCs. 

The main features of ARPES data are determined by 
the dynamical structure of the CMs of Eq. ([T]). In Fig. 
dj we report the spectra along the cuts A-F of Fig. [21 
for x = 0.15 and T = 40 K. We track the QP disper- 
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FIG. 2. (Color online) Fermi surface of LSCO at x = 0.15 
(black solid line). The dashed (red online) and solid (blue 
online) lines mark the C and S hot lines, respectively. The 
shaded circles mark the loci of the waterfalls (from Ref. [l(|). 
The spectra in Figs. Ill3l are calculated along the cuts A-F. 



sions (dots) as the maxima of the momentum distribution 
curves (MDCs) . The condition of quasi criticality (small 
ni\) and the Unite values of the characteristic wave vec- 
tors Qa render the scattering most relevant when the QP 
momenta k satisfy the hot line condition ek = £k±Q A , 
identifying the points at the same energy on the QP 
bands, connected by the characteristic wave vectors. We 
report in Fig. [2] the hot lines for C [(red) dashed lines] 
and S [(blue) solid lines] CMs. These lines intersect the 
Fermi surface at the so-called hot spots. Away from the 
hot lines, the scattering in Eq. [3]is not dominated by the 
nearly critical (low-energy) CM spectrum, and is rather 
mediated by the whole dynamical range of the order of 
V QA . 27j . This determines the energy scale of the kinks 
appearing at low binding energy (< 70meV) in Fig. [1] 
The strong scattering near the hot lines, reminiscent of 
the Bragg scattering occurring when some ordering takes 
place at specific Qa, gives instead rise to the high-energy 
waterfall features in Fig. Q] 33]. 

We obtain waterfalls that compare fairly well with 
the experiments 13, 15, 16 1, although our perturbative 



scheme underestimates their binding energy and broad- 
ening. In particular, our hot lines reproduce the loci of 
the BZ where the waterfalls are observed 16| (shaded 
circles in Fig. [2|). The waterfalls along the cuts A-C (at 
binding energies w 600 meV, « 300 meV, and w 250 meV, 
respectively) correspond to the nearly cross-shaped ac- 
cumulation of the loci well inside the BZ in Fig. [5J In 
our scheme, these are due to C incipient order, which 
also produces additional waterfalls along a square con- 
tour surrounding the T point of the BZ. These are vis- 
ible in panels B and C of Fig. [TJ at approximately 
700 — 800 meV. Their presence cannot be ascertained in 
Ref. (l6j , where the data at higher binding energy are not 
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FIG. 3. (Color online) Same as in Fig. [T] but on a narrower 
energy range. Although the structure of the kinks is due to 
both S and C CMs, the arrows mark the kinks due to (mostly) 
C or S scattering. 



reported. A reanalysis of the data is required to check 
for the presence of these additional waterfalls. 

On the other, hand both the C and S scattering are re- 
sponsible for the dense occurrence of waterfalls near the 
M points (cuts D-F in Fig. [2]). However, as it is clear 
from panels D-F in Fig. [TJ the waterfalls are shifted to 
lower binding energy in this region of the BZ and merge 
with the kinks. Morover, approaching the hot spots, the 
waterfall evolves into a rounding of the QP dispersions, 
with a spectral intensity vanishing as y/uj 0, 0, l35j . 
This rounding is reminiscent of the additional low-energy 



37J, but not in LSCO, pos- 



kinks observed in BSCCO 36, 
sibly due to a lower resolution. 

To follow the evolution of the kinks, we show in Fig. 
|3] the low-energy spectra, along the same cuts of Fig. [U 
In the nodal (TX) direction (cut A), the waterfall is well 
separated from the low-energy kink, which is more closely 
inspected in Fig. @] (a). Here we report the kinked dis- 
persion separately due to C [(red) squares] and S [(green) 
diamonds] CMs, and to the combination of both [(blue) 
circles]. The binding energy of the kink [(blu) arrow in 
Fig. S] (a)] is clearly set by the C CM [(red) arrow], at 
~ 60 — 70 meV, in good agreement with the experimental 
disp ersion for LSCO at the same doping and temperature 
|38| . The markedly propagating character of the C CM 
(phonon-like away from the hot lines 27, 28, 35|) makes 
its contribution to the kink rather sharp. 

We emphasize that the characteristic energy of the C 
CM is extracted from Raman experiments (Tab. QJ and 
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(see Fig. [3]) at energies depending on the CM dispersion 
and on the position of the cut with respect to the hot 
lines. The kink due to the C CM moves to lower ener- 
gies when the cut moves to the region where the C hot 
line approaches the Fermi surface (cuts B and C). On the 
other hand, the S hot lines track rather closely the Fermi 
surface, giving rise to the above mentioned rounding of 
the dispersion at very low energy (< 20meV). Along the 
cuts D-F, the C and S hot lines have intricate structures, 
which makes it difficult to distinguish the role of the two 
CMs in determining the mixed kink-waterfall structures. 
Nonetheless, by switching on and off the C and S cou- 
plings, we can state that the S CM plays the major role 
in this region of the BZ. 



FIG. 4. (Color online) (a) QP dispersion along the FX di- 
rection for LSCO at x = 0.15 and T — 40 K in the pres- 
ence of both C and S CMs (solid circles, blue online) and 
for C (squares, red online) or S (triangles, green online) 
only, (b) Doping evolution of the kink. Both C and S CMs 
are considered, and T = 40 K. The dashed straight lines 
marks the low-energy dispersion. Inset: doping evolution of 
A = (vhe/vf) — 1 (see text). 



is not adjusted here by introducing additional modes at 
suitably chosen phonon frequencies. On the other hand, 
the more diffusive S CM does not fix an energy scale 
and rather renormalizes the bare QP dispersion over a 
broader energy range, affecting the QP velocities far from 
the kink and making the kink more pronounced. Both C 
and S CMs must be simultaneously taken into account 
in order to quantitatively explain the kink. The analysis 
of Raman spectra shows that the interaction mechanism 
switches from S to C CM with increasing x 



17J. This 



characterizes the doping evolution of the QP dispersion 
in the range x = 0.15 — 0.26 38]. In Fig. @] (b), one sees 
that the S-vs-C switching produces no appreciable effect 
on the low-energy dispersion, which is determined coop- 
eratively by the two CMs, so that the slope remains quite 
constant in the doping considered range [dashed line in 
Fig. HI (b)] . On the other hand, in the high-energy disper- 
sion, where most of the effect of the C CM is exhausted 
and the variation of the slope is controlled by the the 
S CM, the dispersion becomes less steep with increasing 
doping, so that the high-energy QP velocity decreases. 
This is clearly observed looking at the inset in Fig. [3] 
(b), where the parameter A = (vhe/vf) — 1, measur- 
ing the strength of the (mostly S-mediated) interaction, 
is plotted vs. x. Although our analysis was limited to 
x > 0.15, extrapolating to the underdoped region the in- 
creasing strength of the S-mediated scattering (see Ref. 
17 1), we can also account for the observed |38[ increase 
of the parameter A (mainly due to the increase of vhe) 
below x — 0.15. 

Along different cuts in the BZ, the QP dispersion is dif- 
ferently affected by the two CMs, which may yield kinks 



The C-S cooperative behavior might be specific of 
LSCO, where the tendency to charge ordering seems to be 
more pronounced than, e.g., in YBCO, where the kinks 
are more rounded. We also stress that our analysis only 
holds above T c . Below T c the S CM changes and dis- 
plays the peculiar resonance at (ir,ir), which alters the 
shape of the kinks, producing a characteristic s-shaped 
dispersion in the antinodal regions [3, H^] . On the other 
hand, the rather broad and moderately coupled phonons 
should keep their effects (most pronounced around the 
nodal regions) even in the superconducting phase. 

In conclusion, the salient aspects of ARPES experi- 
ments in LSCO are well reproduced by the same two 
(C and S) CMs previously obtained to fit Raman exper- 
iments. This solves for LSCO the long-standing issue 
whether the kinks are due to phonons or spin fluctua- 
tions: we reach the Solomonic conclusion that both play 
a role. By the interplay of the two CMs, we can explain 
the highly non-trivial doping evolution of the low- and 
high-energy QP velocity along the nodal (TX) direction, 
with vp almost doping independent and vhe decreasing 
with increasing doping (along with the suppression of the 
coupling with the S CM). We predict the presence of mul- 
tiple kinks (actually, a kink and a low-energy rounding, 



37J). We 



analogous to those observed in BSCCO [3£ 
also predict additional waterfalls at high binding energy, 
along a square contour around the T point of the BZ. 



Since, our analysis demonstrates the presence of two 
CMs, with characteristic wave vectors representative of 
stripe-like textures, our phenomenological model sub- 
stantiates the presence of a competing C and S quasi- 
ordered phase compatible with fluctuating stripes. The 
assessed relevant role of C and S CMs in LSCO also iden- 
tifies them as candidate mediators of the pairing glue in 
these systems. 
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from "University Research Project" of the "Sapienza" 
University n. C26A115HTN. 
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